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Abstract: 

A new measurement of the proton structure function F2(x, Q 2 ) is reported for mo- 
mentum transfers squared Q 2 between 1.5 GeV 2 and 5000 GeV 2 and for Bjorken 
x between 3 • 10 -5 and 0.32 using data collected by the HERA experiment HI in 
1994. The data represent an increase in statistics by a factor of ten with respect 
to the analysis of the 1993 data. Substantial extension of the kinematic range 
towards low Q 2 and x has been achieved using dedicated data samples and events 
with initial state photon radiation. The structure function is found to increase 
significantly with decreasing x, even in the lowest accessible Q 2 region. The data 
are well described by a Next to Leading Order QCD fit and the gluon density is 
extracted. 
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1 Introduction 



A prime task of the electron-proton collider HERA is the investigation of the structure of 
the proton. Measurements of the inclusive lepton-proton scattering cross section have been 
crucial for the understanding of proton substructure p|. Early electron-proton scattering 
experiments have discovered pointlike proton constituents by observing a scale invariant de- 
pendence of the proton structure function F2(x, Q 2 ) on the four-momentum transfer squared 
Q 2 at Bjorken x > 0.1 and Q 2 values of about 5 GeV 2 . Subsequent neutrino scattering ex- 
periments have established the Quark Parton Model (QPM) as a valid picture of the valence 
and sea quarks as constituents of the proton. The interaction of these partons as mediated by 
gluons is successfully described by Quantum Chromodynamics (QCD) which has been tested 
with high precision in muon-nucleon deep-inelastic scattering (DIS) experiments. Experi- 
ments at HERA extend the previously accessible kinematic range up to very large squared 
momentum transfers, Q 2 > 10 3 GeV 2 , and down to very small values of Bjorken x < 10 . 

The first measurements of the structure function F2(x,Q 2 ) reported at HERA, based on 
data collected in 1992, revealed its strong rise at low x < 10~ 2 with decreasing x |2|, ||]. 
This rise was confirmed with the more precise data of 1993 |, §, based on an order of 
magnitude increase in statistics. Such a behaviour is qualitatively expected in the asymptotic 
limit of Quantum Chromodynamics Q. It is, however, not clear whether the rise of F2 is 
fully described by the linear QCD evolution equations, such as the conventional DGLAP 
evolution in logQ 2 or by the BFKL evolution |8|] in log(l/x), or whether there is a 
significant effect due to non-linear parton recombination |J. Furthermore, it is also unclear 
whether this rise will persist at low values of Q 2 of the order of one GeV 2 . For example, Regge 
inspired models expect F2 to be rather flat as function of x at small Q 2 . The quantitative 
investigation of the quark-gluon interaction dynamics at low x is one of the major challenges at 
HERA. It requires high precision for the F2 measurement and complementary investigations 
of the characteristics of the hadronic final state [[II]] . 

In this paper an analysis is presented of inclusive deep-inelastic scattering data taken by 
the HI collaboration in 1994 with an integrated luminosity of 2.7 pb" 1 , which is an order 
of magnitude larger than in 1993. The incident electronP] energy E e was 27.5 GeV and the 
proton energy E p was 820 GeV. The accessible kinematic range has been extended to the 
very high Q 2 region and the structure function F2 has been investigated at a new level of 
precision. To reach lower Q 2 values and correspondingly lower x values, special samples were 
analysed of events with shifted interaction vertex, and of events with tagged initial state 
photon radiation. 

This paper is organized as follows. After a short introduction to the kinematics of inclusive 
ep scattering (section 2), the HI apparatus is briefly sketched (section 3). The different 1994 
data samples, the luminosity determination and the Monte Carlo simulation are described in 
section 4. Next the event selection including the background rejection (section 5) is discussed 
for the different data samples used. Section 6 describes the F2 analyses. In section 7 the 
results are discussed. A phenomenological analysis of F2 is performed and the data are 
compared to recent model calculations at low Q 2 . The data are also studied in the framework 
of perturbative QCD and the gluon distribution is extracted. The paper is summarized in 
section 8. 

X HERA operated with e~p collisions in 1992, 1993 and the start of 1994, and e + p collisions for the major 
part of 1994. In this paper the incident and scattered lepton will always be referred to as an "electron". 
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2 Kinematics 



The structure function F2(x, Q 2 ) is derived from the inclusive electron-proton scattering cross 
section. It depends on the squared four-momentum transfer Q 2 and the scaling variable x. 
These variables are related to the inelasticity parameter y and to the total squared centre 
of mass energy of the collision s since Q 2 = xys with s = AE e E p . A salient feature of the 
HERA collider experiments is the possibility of measuring not only the scattered electron 
but also the complete hadronic final state, apart from losses near the beam pipe. This means 
that the kinematic variables x, y and Q 2 can be determined with complementary methods 
which are sensitive to different systematic effects. These methods were exploited and detailed 
already in Q which describes the analysis of the 1993 data. An appropriate combination of 
the results ensures maximum coverage of the available kinematic range. 

The methods used in the analysis of the 1994 data are the so called "E" (electron) method 
using only the information of the scattered electron and the so called "£" method calculating 
the kinematics based on both the scattered electron and the hadronic final state measure- 
ments [11 1 . The E method, which is independent of the hadronic final state, apart from the 



requirement that the interaction vertex is reconstructed using the final state hadrons, has 
at large y the best resolution in x and Q 2 but needs sizeable radiative corrections. At low 
y the E method is not applied due to the degradation of the y e resolution as 1/y. The £ 
method, which has small radiative corrections, relies mostly on the hadronic measurement 
which has still an acceptable resolution at low y values and can be used from very low to 
large y values. The E and £ results were compared in order to control the calculation of the 
systematic errors. The basic formulae for Q 2 and y for the E method are: 

1 E 'e • 2 2 An/ip 2 #e E 2 SUl 2 6 e 

Ve = 1 - sm — Q e = 4E e E e cos — = — (1) 

bj e I 1 1 — y e 

where E' e and 9 e are the energy and polar angle of the scattered electron. The polar angle 
6 e is defined with respect to the proton beam or z direction, termed "forward" region. The 
formulae for the £ method are 

VS ~ £ + ^(l-COS0 e ) ^ " l-y S 

with 

Z = J2(E h - Pz , h ). (3) 
h 

Here E^ and p z h are the energy and longitudinal momentum component of a particle h, 
the summation is over all hadronic final state particles and the masses are neglected. The 
denominator of is equal to 2E e but measured with all secondary particles. Thus 

ys = ^— 4) 

1 + Vh~ Ve 

with the standard definition 



Vh = — ■ (5) 
The variable x is calculated as x = Q 2 /ys. 



2E e 
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3 The HI Detector 



The HI detector [12] is a nearly hermetic multi-purpose apparatus built to investigate the 
inelastic high-energy interactions of electrons and protons at HERA. The structure function 
measurement relies essentially on the inner tracking chamber system and on the backward 
electromagnetic and the liquid argon calorimeters which will be described here briefly. 

The tracking system includes the central tracking chambers, the forward tracker modules 
and a backward proportional chamber. These chambers are placed around the beam pipe 
at z positions between —1.5 and 2.5 m. A superconducting solenoid surrounding both the 
tracking system and the liquid argon calorimeter provides a uniform magnetic field of 1.15 T. 

The central jet chamber (CJC) consists of two concentric drift chambers covering a polar 
angle range from 15° to 165°. Tracks crossing the CJC are measured with a transverse 
momentum resolution of 5pt/pt < 0.01-jjy/GeV. The CJC is supplemented by two cylindrical 
drift chambers at radii of 18 and 47 cm, respectively, to improve the determination of the z 
coordinate of the tracks. A proportional chamber is attached to each of the z drift chambers 
for triggering. 

A tracking chamber system made of three identical modules measures hadrons emitted in 
the forward direction (7° to 20°). The forward tracker (FT) is used to determine the vertex 
for the events which leave no track in the CJC. This allows an extension of the analysis to 
larger x values. 

In the backward region, attached to the backward electromagnetic calorimeter (BEMC), a 
four plane multiwire proportional chamber (BPC) was located with a polar angle acceptance 
of 151° to 174.5°. The BPC provides a space point for charged particles entering the BEMC 
which is used for low Q 2 < 120 GeV 2 events to identify electrons and to measure 9 e . The 
spatial resolution for reconstructed BPC hits is about 1.5 mm in the plane perpendicular to 
the beam axis. 

The backward electromagnetic calorimeter [O ] which detects the scattered electron at 
low Q 2 is made of 88 lead/scintillator stacks with a size of 16 x 16 cm 2 and a depth of 22 
radiation lengths corresponding to about one interaction length. Around the beampipe the 
stacks are of triangular shape. The angular coverage of the BEMC is 155° < e < 176°. A 
1.5 cm spatial resolution of the lateral shower position is achieved using four photodiodes 
which detect the wavelength shifted light from each of the scintillator stacks. A scintillator 
hodoscope (TOF) situated behind the BEMC is used to veto proton-induced background 
events based on their early time of arrival compared with nominal ep collisions. 

Hadronic final state energies and the scattered electron at high Q 2 (Q 2 > 120 GeV 2 ) are 
measured in the liquid argon (LAR) calorimeter |l4| which covers an angular region between 
3° and 155°. The calorimeter consists of an electromagnetic section with lead absorber 
plates and a hadronic section with stainless steel absorber plates. Both sections are highly 
segmented in the transverse and longitudinal directions with about 44000 cells in total. The 
electromagnetic part has a depth between 20 and 30 radiation lengths. The total depth of 
both calorimeters varies between 4.5 and 8 interaction lengths. 

The luminosity was determined from the measured cross section of the Bethe Heitler 
(BH) reaction ep — > epj. The final state electron and photon can be detected in calorimeters 
(electron and photon "taggers") close to the beam pipe but at large distances from the main 
detector (at z = —33 m and z = —103 m). 
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4 Data and Monte Carlo Samples 



4.1 Data Samples 

Several data samples have been analysed in order to cover maximally the kinematic plane. 
The distribution of the events is shown in Fig. ||. The majority of the events are produced with 
the interaction vertex centered around zero in z, called the "nominal vertex" sample (shown as 
regions C and D in Fig. |lj). Throughout this paper, the low (high) Q 2 sample refers to events 
in which the scattered electron has been detected in the BEMC (L Ar calorimeter) . To reduce 
the systematic errors of the Fi measurement, a strict data selection was performed based on 
the behaviour of the main detector components. This behaviour was required to be optimal 
for the low Q 2 analysis of the nominal vertex sample which allows the highest precision to 
be reached. The remaining integrated luminosity for the low Q 2 sample is 2.2 pb _1 , the one 
for the high Q 2 sample is 2.7 pb _1 . The number of accepted events per unit luminosity was 
checked to be constant within statistical errors during the data taking period. 

In order to study the behaviour of F2 at small Q 2 several means were used to extend the 
acceptance to this kinematic region using special event samples. For DIS events at very low 
Q 2 the electron is scattered through a large angle 9 e . For 9 e values greater than 173° and the 
interaction vertex at its nominal position at z = +3 cm the electron hits the inner edge of 
the BEMC calorimeter or remains undetected near the beam pipe. The acceptance extension 
in the backward region was realized as follows: 

• During good accelerator background conditions, the innermost parts of the backward 
electromagnetic calorimeter (BEMC) around the beampipe were included in the trigger 
for part of the time. Since these detector elements are of triangular shape, these data 
will be referred to as the "open triangle" data sample. An integrated luminosity of 
0.27 pb _1 was accumulated. The kinematic region covered by this sample is shown as 
region (B) in Fig. |]. 

• As in 1993 Q, the interaction point was shifted in the forward direction to an average 
position of z = +67 cm which permits measurements up to 9 e ~ 176.5°. This sample 
of 58 nb _1 of data is referred to as the "shifted vertex" data sample to distinguish it 
from the data with a nominal event vertex. It covers region (A) in Fig. |l]. 

• The low Q 2 region was also accessed by analyzing events from the so called early 
proton satellite bunch colliding with an electron bunch at z ~ +68 cm. The kinematic 
region covered by this sample is similar to that of the shifted vertex data sample. The 
"satellite" data sample amounts to ~ 3% of the total data corresponding to a total 
"luminosity" of 68 nb _1 selected over the whole run period. 

• Finally, a sample of deep-inelastic radiative events was extracted with a hard photon 
emitted collinear with the incident electron. These events have a reduced incident 
electron beam energy which allows access to very low Q 2 values with the present detector 
setup. Since only about 2% of the DIS events are tagged as radiative events, the nominal 
vertex sample had to be used for this study. Subsequently the tagged radiative events 
are referred to as "the radiative event sample" and the bulk of the data are sometimes 
called "non-radiative" in contrast. 



7 



HI 1994 




10 100 1000 

Q 2 [GeV 2 ] 



Figure 1: Distribution of the event sample in the (x,Q 2 ) plane. The 4 visible regions 
(A,B,C,D) correspond to A) events recorded during a period in which the interaction region 
was shifted with respect to the nominal position allowing access to larger 9 e ; B) events 
from the nominal vertex position taken in a period in which the innermost BEMC stacks 
of triangular shape were included in the trigger ("opened triangles", see text) or C) not 
included; D) high Q 2 events with the scattered electron detected in the LAr calorimeter. 
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4.2 Luminosity Determination 



The most precise method of determining the luminosity from the reaction ep — > ep^y is 
based on the measurement of the energy spectrum of hard photons (JSL > 10 GeV) as 
explained in [15] for the 1993 data. The main uncertainties of the measurement of the 



integrated luminosity for the 1994 nominal vertex data are: the photon tagger absolute 
energy scale (0.9%), the trigger efficiency of the luminosity system (0.3%), the precision of the 
electron gas background subtraction (0.4%), the photon-tagger acceptance (0.5%), multiple 
photon overlaps (0.4%), the precision of integration resulting from the 10 sec interval between 
consecutive luminosity measurements (0.5%) and the correction for satellite bunches (0.5%). 
Major improvements with respect to 1993 data include the trigger efficiency, the satellite 
bunch correction and the precision of the energy scale in the photon tagger. The precision of 
the luminosity measurement for the nominal vertex data sample is 1.5% which represents an 
improvement of a factor of 3 with respect to the 1993 data analysis. For the shifted vertex 
data sample the luminosity uncertainty is 3.9%. 

The results of this measurement were checked for consistency with a sample of Bethe 
Heitler events in which both the electron and photon are detected simultaneously, and with 
QED Compton events. Both these analyses are subject to different systematics, compared 
with the hard photon method, allowing a cross check of the luminosity with a precision of up 
to 6%. 

The integrated luminosity of the satellite data sample was obtained from the measured 
integrated luminosity for the shifted vertex data multiplied by the efficiency corrected event 
ratio in a kinematic region common to both data sets. The precision of that luminosity 
determination was estimated to be 7.1%. 



4.3 Monte Carlo Simulation 



More than one million Monte Carlo DIS events were generated using the DJANGO 16] pro- 
gram. The Monte Carlo event statistics correspond to an integrated luminosity of approxi- 
mately 18 pb" 1 . The DJANGO program is based on HERACLES fl7[ | for the electroweak 
interaction and on the LEPTO program [|l^] to simulate the hadronic final state. HERACLES 
includes first order radiative corrections, the simulation of real Bremsstrahlung photons and 
the longitudinal structure function. The acceptance corrections were performed using the 



GRV parametrization [19] which describes rather well the HERA F% results based on the 
1993 data. LEPTO uses the colour dipole model (CDM) as implemented in ARIADNE @ 
which is in good agreement with data on the energy flow and other characteristics of the 
final state as measured by HI |21] and ZEUS [22|. For the determination of systematic errors 



connected with the topology of the hadronic final state, the HERWIG model 23] was used. 



Photoproduction background was simulated based on the PHOJET p4|] , PYTHIA ]2J 
and RAYVDM f26f| generators for jp interactions. With these models large samples of photo- 
production events were generated which contained all classes of events (elastic, soft hadronic 
collisions, hard scattering processes and heavy flavour production). 

It was found that about 10% of the DIS data at HERA consists of events with a large 
gap in pseudo-rapidity around the proton remnant direction [p7j] . These events were found to 



be compatible with diffractive exchange and are well described by the model RAPGAP [28] 
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as deep-inelastic scattering on a colourless object -termed a pomeron- emitted from the 
proton. The hadronic final state of these events is also well described by RAPGAP which 
includes ARIADNE for QCD effects. The RAPGAP Monte Carlo simulation was used to 
check the effect of the large rapidity gap events on the vertex reconstruction efficiency which 
depends mostly on the final state topology of the events. Differences between rapidity gap 
and "standard" DIS events of up to 2% were found at large y > 0.4 and smaller at low y, 
and were included in the systematic error of i^. 

For the events generated with the models described above the detector response was 



simulated in detail [12] using a program based on GEANT [29|. The simulated Monte Carlo 



events were subjected to the same reconstruction and analysis chain as the real data. 



5 Event Selection 



The low Q 2 DIS events in the backward region were triggered by an energy cluster in the 
BEMC (E' e > 4 GeV) which was not vetoed by the TOF. The high Q 2 events were triggered 
by requiring an electromagnetic energy cluster in the LAr calorimeter (E' e > 8 GeV) . A trigger 
of lower energy threshold {E' e > 6 GeV) also accepted the event if there was simultaneously 
a tracking trigger. In the region of the final F2 data presented below the trigger efficiency, 
which has been determined from the data, is about 80% for E' e ~ 8 GeV, and becomes larger 
than 99% for E' e > 10 GeV. 



5.1 Selection of Deep-Inelastic Scattering Events 

Deep-inelastic scattering events in HI are identified by the detection of the scattered electron 
in the BEMC or LAr calorimeter and the presence of a reconstructed interaction vertex. The 
electron identification cuts, fiducial volume and vertex requirement are detailed in Table 1. 

These selection criteria follow closely those of the 1993 data analysis Q. For the low Q 2 
nominal vertex sample (Q 2 < 120 GeV 2 ) an additional cut tbpc < 64 cm is applied, where 
r BPC is the radial distance of the electron hit in the BPC to the beam axis. This cut prevents 
the electron from entering the transition region between the BEMC and the LAr calorimeter 
where the energy corrections are large and depend strongly on the impact point. For the 
same reason, the high Q 2 events (Q 2 > 120 GeV 2 ) are accepted only if the electron cluster 
is fully contained in the LAr calorimeter. Despite these conditions, the measurement could 
also be performed for intermediate Q 2 (Q 2 ~ 120 GeV 2 ) due to the ±30 cm spread of the 
event vertex position around its nominal position. 

The scattered electron is identified with the electromagnetic cluster of maximum energy 
which satisfies the estimator cuts of Table 1 . The electron identification efficiency, determined 
from Monte Carlo simulation studies, is better than 97% except at Q 2 < 6.5 GeV 2 where it 
falls to 94% at the lowest x values. 

At low Q 2 the main sources of non-ep background are due to proton beam interactions 
with residual gas and beam line elements upstream of the HI detector. At high Q 2 the main 
background is due to cosmic ray events and muons travelling off axis parallel to the proton 
beam. An efficient reduction of these background contributions is provided by the minimum 
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low ^smrtecl vtxj 


low Q 2 (nominal vtx) 


high Q 2 




H J TY1 Dl rl (H 

Jjj, Zj IIlcLIlUU. 


H » m of l~i y~\ f~\ 
JJJ, Zj IIlcLIlUU. 


E method 


£ method 


e e /° 


< 176 


< 173 


< 150 


< 153 


E'JGeV 


> 11 


> 11 


> 11 


> 11 


^vertex /CIH 


67± 30 


5± 30 


5± 30 


5± 30 


electron identif. 


ei < 5 cm 


t\ < 5 cm 


e 3 > 50% 


e 3 > 65% 


electron identif. 


£2 < 5 cm 


£2 < 5 cm 


e 4 > 3% 


es < 30 mrad 



Table 1: Summary of event selection criteria for the shifted and the nominal vertex (vtx) 
data at low and high Q 2 . The approximate event numbers are 10000, 220000 and 9000 
events respectively. For the open triangle data sample the 9 e cut is 174°. For the electron 
identification several estimators were used: e±: electron cluster radius; smallest distance 
from the closest hit in the BPC to the centroid of the electron cluster; e 3 : fraction of the 
electron energy deposited in the 4 most energetic cells of the cluster; e^: fraction of the 
electron energy deposited in the first 3 radiation lengths of the calorimeter; e^: angle between 
the line connecting the vertex to the centroid of the electron cluster and the associated track. 



energy and the vertex requirements discussed above. The number of residual beam-induced 
background events was estimated from non-colliding bunch studies, and the number of cosmic 
events from scanning. Both together represent less than 1% of the number of selected events 
in any (x, Q 2 ) bin. 

The only significant background to DIS from ep interactions is due to photoproduction 
events where the scattered electron escapes the detector along the beam pipe but in which an 
energy cluster from the hadronic final state fakes a scattered electron. About 10% of these 
events are identified as photoproduction background if the scattered electron is found in the 
electron tagger. Photoproduction events were simulated to estimate this background. The 
photoproduction background was subtracted statistically bin by bin. Only 12 bins, out of a 
total of 193 (x,Q 2 ) bins, have a contamination larger than 3%. This contamination never 
exceeds 15% in any bin. 

Figure ||a shows the distribution of the angle of the scattered electron for the shifted 
vertex data compared to the Monte Carlo simulation weighted with the measured structure 
function (see section 6). The Monte Carlo simulation is normalized to the luminosity and 
agrees well with the data illustrating the level of residual background in the low Q 2 sample. In 
Fig. |2]b the normalized energy spectrum in the electron tagger is shown which is reproduced 
by the background photoproduction event simulation. 

Figure ||a shows the distribution of the energy of the scattered electron for the high 
statistics nominal vertex data. The simulation gives an excellent description of the data from 
the low energy up to the so called kinematic peak region, i.e. the region around the value of 
the incident electron beam energy. This agreement was achieved after a spatially dependent 
calibration of the data and Monte Carlo response [ 30 1 using the double-angle method [31]. 



The small remaining contribution of the photoproduction background is also shown. In 
Fig. |3]b the fractions of y^ originating from tracks, BEMC and LAr calorimeter are given as 
a function of log 10 yn- In this analysis the yh variable is determined by using a combination 
of central tracks and calorimeter cells An isolation criterion is used to avoid counting 
the energy of the LAr cells originating from a track already used in y^. For y <0.15, i.e. in 
the region where the E method will be used for the F2 result, the measurement is dominated 
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Figure 2: Shifted vertex data: experimental and Monte Carlo distributions of a) the polar 
angle of the scattered electron and b ) the energy of the scattered electron in photoproduction 
background events detected in the electron tagger. 




log (y h ) 



Figure 3: Nominal vertex data with the scattered electron in the BEMC (Q 2 < 120 GeV 2 ): 
experimental and Monte Carlo distributions a) of the scattered electron energy and b) of the 
fraction of contributed by the tracks, the LAr calorimeter and the BEMC 



12 
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Figure 4: High Q 2 data: experimental and Monte Carlo distributions a) of the energy of the 
scattered electron detected in the LAr calorimeter and b) of the ratio y^/y e , for y e > 0.05. 

by the track reconstruction and the LAr measurement (Fig. ^b). At larger y, the BEMC 
contribution plays an increasing role due to the low energy particles which accumulate in the 
backward direction. The DIS Monte Carlo simulation describes well these fractions in the 
complete kinematic range. 

Figure |I|a shows the distribution of the energy of the scattered electron detected in the LAr 
calorimeter. It is well described by the Monte Carlo simulation. A detailed calibration was 
carried out by comparing events from the kinematic peak at low y (< 0.1) with simulation, 
including corrections for the energy lost due to the dead material between the wheels which 



make up the LAr calorimeter 32]. This procedure has been cross checked with the double- 
angle method. 

Figure ||]b shows the ratio ys/y e i n the high Q 2 sample compared to the Monte Carlo 
expectation. The resolution of this ratio, which is calculated for y e > 0.05, and thus of 
is better than 13% in this kinematic region. The "tail" visible at values below 0.7 is due to 
radiative events, and is well described by the Monte Carlo simulation. 



5.2 Selection of Deep-Inelastic Radiative Events 

A sample of deep-inelastic events with an energetic photon {E~ > 4 GeV) emitted collinear 
with the incident electron was selected. These radiative events can be interpreted as deep- 
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Figure 5: Radiative events: experimental and Monte Carlo distributions of A (eq. 6) with 

a) the energy detected in the electron tagger (E etag ) bigger than 2 GeV; b) with E etag < 2 
GeV and c) distribution of the photon energy detected in the photon tagger. The analysis 
cut in b) indicates the region of A > 0.5 excluded from the analysis. The full solid line in 

b) and c) represents the sum of all three contributions in the Monte Carlo: DIS initial state 
radiation events (ISR MC), DIS events with a BH overlap (DIS + BH) and photoproduction 
events with a BH overlap (jp + BH). 



inelastic scattering events with a reduced ("true") incident energy E± = E e — E~ which can 
be reconstructed due to the additional detection of the radiated photon in the small angle 
photon tagger of the luminosity system. When using the E method, the kinematic variables 
yt and Q 2 are obtained by replacing in eq.[l] the nominal beam energy by the reduced energy 
Ef. Note that Q 2 ^ and are unchanged by the E e — ► Et transformation while xs is affected. 

A first experimental study of this process at HERA has been published ||, H by the HI 
collaboration using 1993 data, which where however too limited in statistics for a quantitative 
study of the proton structure. The larger integrated luminosity of the 1994 data permits a 



significant F2 measurement for Q 2 values down to 1.5 GeV 2 . The ZEUS collaboration [36] 
recently published results on F2 using this method. 

A summary of the selection criteria of the final sample of about 8200 events is given in 
Table ||| |35|]. The event selection for radiative events is similar to the one for low Q 2 non- 
radiative events, apart from the additional requirement of a detected photon with at least 
4 GeV in the small angle photon tagger of the luminosity system. This requirement also 
reduces the photoproduction background. Therefore the minimum scattered electron energy 
can be lowered to 8 GeV. 

The selected sample contains both radiative DIS events and pile-up events due to overlaps 
of DIS and "fp events with Bethe Heitler events in a time window of ±5 ns. The pile-up events 
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are partly removed from the sample by requiring the energy in the electron tagger, E e t ag , to 
be less than 2 GeV, but the majority of them remains. 

The background can be controlled through the redundancy of the true electron beam 
energy measurement Et- For radiative DIS events we expect measurements of the quantity 

A = [S 7 - E e (y e - y h )]/E 1 (6) 

to be concentrated around zero while for pile-up DIS events a concentration around one is 
expected. Here y e and yh are calculated according to eqs. [j] and |B[ The distribution of A 
is shown in Fig. ^ for a sample of events with a) E etag > 2 GeV (dominantly ep collisions 
with BH overlap events) and b) E e t ag < 2 GeV. The data are compared with Monte Carlo 
simulation. The pile-up sample in Fig. |5]a shows a clear peak for A = 1, and is well described 
by the sum of simulated DIS and ^yp distributions with overlap of BH events. Fig. |B|b shows 
a peak for A = 1 from residual pile-up events for which the electron from the BH event was 
not detected, and a peak around A = from genuine radiative events. Radiative events are 
selected in this analysis by requiring A < 0.5. The background of pile-up events as estimated 
by the Monte Carlo simulation studies is subtracted statistically. The remaining background 
from overlap 7p and DIS events estimated from Monte Carlo studies amounts to 8%, with at 
most 15% in an x, Q 2 bin. In Fig. ||c the photon energy spectrum as measured in the photon 
tagger is shown for the selected sample and compared with simulated signal and background 
distributions. There is a good agreement between data and Monte Carlo simulation. 





low Q 2 (radiative events) 




E, £ method 


e /° 


< 174 


E'JGeV 


> 8 


^vertex / Cm 


5± 35 


electron identif. 


e\ < 5 cm 


electron identif. 


62 < 4 cm 


£ 7 /GeV 


> 4 


E e t ag /GeV 


< 2 


A 


< 0.5 



Table 2: Summary of event selection criteria for the radiative event sample. For the electron 
identification two estimators were used: e±: electron cluster radius and e%: smallest distance 
from the closest track to the centroid of the electron cluster. The variable A is dehned in 
eq. |. 



6 Structure Function Measurement 



The structure function ^(x, Q 2 ) was derived from the one-photon exchange cross section 



d 2 a 



27TQ 2 



2y + 



y 



;)E2(x,Q 2 ). 



(7) 



dxdQ 2 Q A x v " ' 1 + R 

The structure function ratio R = F2/2xF\ — 1 has not yet been measured at HERA. It was 
calculated using the QCD relation [37] with the NLO strong coupling constant [38] and the 
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GRV structure function parametrization. Note that a 20% error on R corresponds to about 
2% uncertainty on F2 at y = 0.6 for R of about 0.6. The R values are quoted in Tables 
and H| no extra effect of the R uncertainty on F2 was considered. 

Compared to the previous HI analysis the F2 measurement has been extended to 
lower and higher Q 2 (from 4.5 — 1600 GeV 2 to 1.5 — 5000 GeV 2 ), and to lower and higher 
x (from 1.8 • 10 -4 < x < 0.13 to 3 • 10~ 5 < x < 0.32). The determination of the structure 
function requires the measured event numbers to be converted to the bin averaged cross 
section based on the Monte Carlo acceptance calculation. The binning in x was governed 
by the detector resolution and could be chosen to be rather fine since the E and £ methods 
were used in the optimum range at low and high x, respectively. The x resolution is better 
than 20%. The Q 2 resolution is about 5% and the number of bins in Q 2 was adapted to 
statistics. All detector efficiencies were determined from the data utilizing the redundancy 
of the apparatus. Apart from very small extra corrections, all efficiencies were correctly 
reproduced by the Monte Carlo simulation. The bin averaged cross section was corrected for 



higher order QED radiative contributions using the program HECTOR [39]. Effects due to 
Z boson exchange at present values of Q 2 and y are smaller than 3% and were treated as 
part of the radiative corrections. 

Different data sets are available which, for a given (Q 2 ,x) interval, use different parts of 
the detectors. Thus many cross checks could be made in kinematic regions of overlap for the 
two kinematic reconstruction methods and these gave very satisfactory results. In this paper 
results are presented from the radiative F2 analysis (1.5 < Q 2 < 3.5 GeV 2 ), from the shifted 
vertex analysis (1.5 < Q 2 < 2.5 GeV 2 ), from a combination of the shifted vertex and the 
satellite bunch analysis (3.5 < Q 2 < 6.5 GeV 2 ), from the open triangle analysis (Q 2 = 8.5 
GeV 2 ) and from the nominal high statistics sample when the scattered electron is detected 
in the BEMC (12 < Q 2 < 120 GeV 2 ) or in the LAr calorimeter (120 < Q 2 < 5000 GeV 2 ). 

Compared with the analysis of the 1993 data, many uncertainties have been reduced. The 
systematic errors are due to the following sources: 



• The uncertainty in the electron energy scale which is 1% in the BEMC, and 3% in the 
LAr calorimeter. Since the y e resolution varies as 1/y with the energy resolution even 
a 1% error on 5E/E can lead to 10% errors on F2 at low y in the E method. 

• The uncertainty in the hadronic energy scale: the detailed study of yh/Ue and of 
PT,h/PT,e (pt is the momentum transverse to the beam axis) allowed the assignment of 
a 4% error on the hadronic energy deposited in the LAr calorimeter, a 15% error on the 
same quantity in the BEMC, and a 3% error on the yh fraction carried by the tracks. 
These errors take into account the intrinsic energy scale uncertainty of each detector 
and the uncertainty of the sharing of the total hadronic final state energy between these 
three subdetectors. These numbers also include uncertainties due to the treatment of 
the electronic noise in the LAr calorimeter and the BEMC. 

• An uncertainty of up to 1 mrad for the electron polar angle which leads to an error on 
F 2 of 8% at low Q 2 . 

• Apart from the electron identification, all efficiencies were determined from the data 
and compared with Monte Carlo simulation. The agreement between the experimental 
and the simulated values for the individual efficiencies was found to be better than 
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2%. An overall error of 2% was assigned due to the imperfect description of the var- 
ious efficiencies. A larger error was added to account for the variation of the vertex 
reconstruction efficiency at large x (up to 8%) where jets get closer to the beam pipe 
in the forward direction, and at small x or large (up to 4%) where HI had no further 
tracking device besides the BPC to monitor the vertex efficiency. 

• Uncertainties in the hadronic corrections, the cross section extrapolation towards Q 2 = 
and higher order corrections, which give an error of up to 2% in the radiative correc- 
tion. The accuracy was cross checked by comparing the HECTOR calculation with the 
HERACLES Monte Carlo simulation results. The agreement to the few percent level 
between the structure function results obtained with the E and the £ methods is an 
additional cross check for the control of the radiative corrections. 

• The structure function dependence of the acceptance which was kept below 1% by 
performing a two step iterative analysis. The uncertainty in the simulation of the 
hadronic final state reflects most prominently in the efficiency for the requirement of 
an interaction vertex from tracks. A comparison of the different models (sect 4.3) for 
the hadronic final state was used to assign an additional 3% systematic error entering 
in all analyses at low x through the vertex efficiency. 

• Based on the control data sample of electron tagged 7p events the uncertainty due 
to photoproduction background could be estimated to be smaller than 30% of the 
correction applied. This is equivalent to at most a 5% systematic error in the highest 
y bins at lower Q 2 only. 

• Statistical errors in the Monte Carlo acceptance and efficiency calculations were com- 
puted and added quadratically to the systematic error. 

• For the analysis of radiative events an additional 1.5% uncertainty on the photon energy 
measurement in the photon tagger was considered and a 2% systematic error was added 
due to the uncertainty of the photon tagger geometrical acceptance. An uncertainty on 
the trigger efficiency of 6% to 9% was included for the lowest x points. 

Some of the systematic uncertainties affect differently the Fi measurement in the different 
methods. The systematic errors are given in Tables 7 and 8 point by point. However, some 
are strongly correlated over a large kinematic range. These correlations have been considered 
in the fits reported below. The matrix with the many different error contributions is available 
upon request to the HI collaboration. In Fig. ^ the comparison of the measurements using 
the E and using the X method is shown. Both are in good agreement for all Q 2 values. Some 
possible discrepancies between both methods, e.g. at Q 2 = 3.5 GeV 2 , were investigated 
carefully and taken into account when evaluating the systematic error of the final structure 
function values if they could not be resolved. Small deviations are possible though, due to the 
different population of the x, Q 2 plane between the methods of calculating the kinematics. 
The two measurements are combined using the E method for y > 0.15 and the X method 
for y < 0.15. The result is shown in Fig. [7| and given in Tables 7 and 8. The measurements 
obtained from the low Q 2 nominal vertex sample have a typical systematic error of 5%. The 
large statistics allow the measurement of F2 to reach Q 2 values of 5000 GeV 2 , and to achieve 
a few percent statistical precision at Q 2 below 100 GeV 2 . The results are in good agreement 
with the previous HI publication [|]. In particular, the distinct rise of F2 towards low x, 
observed with the 1992 and 1993 |2 , |I| data, is confirmed with higher precision and extends 
now to significantly lower Q 2 values. 
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Figure 6: Measurement of the structure function with the electron (closed circles) and the 
£ method (open squares). The inner error bar is the statistical error. The full error bar 
represents the statistical and systematic errors added in quadrature disregarding the error 
from the luminosity measurement. 
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Figure 7: Measurement of the structure function F2(x,Q 2 ) as a function of x. The closed 
circles are the results of this analysis, the open circles are results taken from the recent 



publication of the NMC l 4C] and the open squares results from BCDMS The inner error 
bar is the statistical error. The full error bar represents the statistical and systematic errors 
added in quadrature and disregarding the luminosity error. Additionally a data point has 
also been measured at Q 2 = 2 GeV 2 (see Table [?J not shown in the Figure). The curves 
represent the NLO QCD fit described in section 7.4, which includes the data for Q 2 > 5 
GeV 2 . The extension of the curves below 5 GeV 2 represents only the backward evolution of 
the fit. 
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7 Discussion of the Results 



7.1 Phenomenological Fits to the Data 



The x and Q 2 behaviour of F 2 can be described by a phenomenological ansatz of the type 

F 2 {x, Q 2 ) = [a-x b + c-x d -{l + e-^x~)- (log Q 2 + / log 2 Q 2 + h/Q 2 )} ■ (1 - x) 9 , (8) 

where Q 2 is given in GeV 2 . This functional form was introduced in detail previously j|. An 
extra 1/Q 2 term has been added in order to get a good description at Q 2 below 5 GeV 2 . 
Note that this term is not a measure of higher twist contributions. For the fit, results from 



a 


b 


c 


d 


e 


f 


g 


h 


3.10 


0.76 


0.124 


-0.188 


-2.91 


-0.043 


3.69 


1.40 



Table 3: Parameters of a phenomenological fit to the proton structure function results 
from this experiment combined with F 2 measurements from the NMC and the BCDMS 
experiments. The parametrization is valid for 1.5 GeV 2 < Q 2 < 5000 GeV 2 , 3-10~ 5 < x < 1 
and Q 2 < x-10 5 GeV 2 . The parameter h is given in GeV 2 . 



HI, NMC [[5]] and BCDMS [41] are used and statistical and systematic errors were added in 
quadrature. The parameter values quoted in Table ^ are close to those obtained with 1993 
data. The fit provides a good description of all data from the experiments with a x 2 /dof of 
1.65 using full errors. For the HI data alone the parametrization gives a X 2 /dof of 1.00. 



In perturbative QCD the rate of growth of F 2 towards low x is expected to increase with 
increasing Q 2 ||. The wide range of momentum transfer covered in this experiment enables 
a study of the Q 2 dependence of the power A characterizing the rise of F 2 oc x~ x at low 
x. For each Q 2 bin and x < 0.1 the exponent A was determined taking into account the 
point-to-point systematic error correlations. The result is given in Table |] and displayed in 
Fig. H A rise of A with log Q 2 is observed in the covered range from about 0.2 to 0.4 between 
1.5 and 800 GeV 2 . 



The structure function F 2 is related to the total cross-section of the virtual photon-proton 
interaction, (Jtot{l*p), via 

(9) 



At low x, W is equal to y/Q 2 / x - The A parameter thus also determines the dependence 
of F 2 on the invariant mass squared W 2 of the virtual photon-proton {^f*p) system. For 
hadronic and real photoproduction total cross sections the value of A has been measured to 
be around 0.08 [42], which is interpreted as the intercept of the so called soft pomeron. For 



virtual photon-proton interactions A is found to be substantially larger, and increases with 
Q 2 . Future analyses of HERA data, which will lead to F 2 measurements at Q 2 below 1 GeV 2 , 
should allow the transition between deep-inelastic scattering and real photoproduction to be 
studied. 
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Figure 8: Variation of the exponent A from fits of the form F2 ~ x A at fixed Q 2 values and 
x < 0.1. 
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H1 1994 




Figure 9: Measurement of the proton structure function F2(x,Q 2 ) in the low Q 2 region 
by HI (closed circles: non-radiative events; closed triangles: radiative events), together with 
results from the ZEUS |36|/ (open squares), E665 f^j (open circles) and NMC pQ/ (open 
triangles) experiments. The Q 2 values of the ZEUS data shown for the bins Q 2 =3.5, 5 and 
6.5 GeV 2 are measurements at 3.0, 4.5 and 6 GeV respectively. Different parametrizations 
for i*2 are compared to the data. The DOLA and CKMT curves are only shown for the upper 
row of Q 2 bins; CTEQ3M and MRSG are shown for the lower row; GRV and MRSA' are 
shown for the full Q 2 range. 
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8.5 


0.228 


0.015 


0.014 


200 


0.372 


0.039 


0.032 


12 


0.238 


0.007 


0.031 


250 


0.360 


0.036 


0.060 


15 


0.261 


0.006 


0.022 


350 


0.270 


0.063 


0.060 


20 


0.268 


0.007 


0.020 


500 


0.460 


0.108 


0.047 


25 


0.286 


0.009 


0.024 


650 


0.391 


0.128 


0.059 


35 


0.331 


0.010 


0.022 


800 


0.503 


0.235 


0.129 



Table 4: The values of the exponent X as a function of Q 2 . 



7.2 Comparison with Models at Low Q 2 

Figures [?] and || clearly demonstrate the rise of F2 with decreasing x. In Fig. |9] the data from 
the eight lowest Q 2 bins are shown^Jand compared with recent data and F2 parametrizations. 
The rise of F2 towards low x is also present in the low Q 2 region. The measurement is in 
good agreement with the data from the ZEUS experiment |36"| and matches well with the 
data from fixed target experiments [43, ^] at higher x values. 

The curves denoted as MRSA', MRSG g§, CTEQ3M g| and GRV g| are parametriza- 
tions based on the conventional QCD evolution equations. These calculations assume a cer- 
tain shape of the x behaviour at a starting Qq value and use the DGLAP Q equations to 
get predictions at different Q 2 values. The MRS and CTEQ distributions assume an x~ x 
behaviour for x — » at starting Qq of a few GeV 2 . Their parameters were determined using 
also the 1993 HERA structure function data. 

The GRV calculation assumes that all parton distributions at very low Qq = 0.34 GeV 2 
have a valence like shape, i.e. vanish for x — > 0. Assuming that the DGLAP equations can 
be used to evolve the parton distributions from this low Qq scale to larger Q 2 values, they 
predicted that the structure function F% should rise towards low x even for low values of 
Q 2 ~ 1 GeV 2 jit]]. The determination of the shape parameters of the distributions at the 
starting scale uses only data from fixed target experiments and not much freedom is left for 
further adjustments in the kinematic range of the HERA data. Small variations are connected 
with changes still possible in the starting Qg and the value of the QCD parameter A. Fig. || 
shows that the GRV distributions describe the data well, indicating that in this kinematic 
regime the sea quark distributions can be produced by QCD dynamics. 

Parametrizations motivated by Regge theory relate the structure function to Reggeon 
exchange phenomena which successfully describe the slow rise of the total cross section with 



2 For the radiative data points (triangles in Fig. the y variable cannot be calculated using Q 2 jsx with the 
nominal s since each bin has a different mean incoming electron energy. The average y values are 0.143, 0.063, 
0.026, 0.010 for Q 2 = 1.5 GeV 2 and 0.199, 0.086, 0.036, 0.015 for Q 2 = 2.5 GeV 2 , starting at the smallest 
given x value. For the single point at Q 2 — 3.5 GeV 2 y = 0.018. 



23 



the centre of mass system energy in hadron-hadron and jp interactions. Using the "bare" 
instead of the "effective" pomeron intercept, the CKMT |^] parametrization rises faster with 
x compared to former DOLA E7j calculations. The CKMT curves were calculated using the 



parameters as given in [46], without QCD evolution in the whole range. 



The predictions for the Regge inspired models DOLA and CKMT lie below the data for 
Q 2 > 2 GeV at low x. The latter were already shown to be significantly below the HI data of 
1993 pq| . The GRV and MRSA' parametrizations give a good description of the data in the 
range shown, with the possible exception of the first Q 2 bin for the latter. The MRSG and 
CTEQ3 distributions, which are not available for the lowest Q 2 values, describe the higher 
Q 2 data well. 



7.3 Double Asymptotic Scaling 

The success of the GRV approach suggests that the observed rise of the structure function 
F2 towards low x is generated by QCD dynamics. This was already observed in 1974 || from 
a study of the behaviour of F% in the limit of large Q 2 and low x. In this asymptotic region 
the QCD evolution determines the shape of F2. Recently Ball and Forte |^9| developed a 
convenient way to test the asymptotic behaviour of F2 using two variables 



a = ^J\og(x /x) ■log(a s (Q )/a s (Q)), P = ]J 



log{x /x) ^ 



log(a s (Q )/a s (Q)) 



where a s {Q) is evaluated at the two loop level [50|. 

The parameters xq and Qq have to be determined experimentally. The parameter Qq 
is optimized by minimizing the \ 2 °f a linear fit of log(R' F F2) versus a (see below) using 
data with Q 2 > 5 GeV 2 . This leads to a value of Qq = 2.5 GeV 2 . The same procedure was 



followed for xq, which showed less sensitivity. The value xq = 0.1, as suggested in |49|, 5jj , 
was found to be a good choice. To visualize the double scaling, it was proposed to rescale F2 
with factors R' F and Rf defined as 

R F (a, p) = 8.1 exp f-2 7 a + u- + - logger) + log(^)) /£p (11) 
V P 2 7 J 



with 
and 



tF = 1 + ((6 + 6) * <*s(Q) ~ fi * «»(0o)) * (p/(2tt * 7)) (12) 



R' F (a,p) = R F exp{2ja). (13) 

Here £1 = (206n//27 + 6b 1 /b )/b , & = 13, b = 11 - 2n f /3, lo = (11 + 2n//27)/6 and 
b\ = 102 — 38nj/3. The number of flavours is nf and 7 = >/ (12/6o)- The function log(R' F F2) 
is then predicted to rise linearly with a. RfF2 is expected to be independent of p and a. 
Note that these expectations are valid only if the gluon distribution, which drives F2 at low 
x via the sea quarks, does not have a too singular behaviour for Q 2 = Q\. 

Fig. [To| a shows RpF2 versus p for the data with Q 2 > 3.5 GeV 2 . The value of A for 



four flavours was chosen to be 263 MeV [52|. The continuity of a s (Q) at the bottom quark 
mass threshold is imposed using the prescription in [j38| l. Approximate scaling is observed for 
Q 2 > 5 GeV 2 and p > 2. At high p the low Q 2 data tend to violate the scaling behaviour 
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Figure 10: The rescaled structure functions a) RfF2 versus p and b) log(R' F F2) versus a 
(see text). Only data with Q 2 > 5 GeV 2 and p > 2 are shown in b). 
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which is seen clearly for the data at 3.5 GeV 2 . In Fig. 10b, the results are shown for p > 2 and 
Q 2 > 5 GeV 2 as a function of a. The data exhibit the expected linear rise of log(i?^i 7 2) with 
a. A linear fit to the data gives a value for the slope of: 2.50±0.02±0.06 (2.57±0.05±0.06) for 
Q 2 < 15 GeV 2 (Q 2 > 35 GeV 2 ) and 4 (5) flavours. The first error is the statistical error and 
the second error is the systematic error taking into account the point-to-point correlations. 
The value expected from QCD is 2.4 (2.5) for 4 (5) flavours. The results are in agreement 
with these predictions. Compared to the result presented in [pi]] , the extraction based on 
the 2-loop formalism used here is in better agreement with QCD expectation. Not included 
in this error is the influence of the uncertainty in the choice of A. Varying A by ±65MeV 
changes the result on the slope by =f0.03. 

One can conclude that the low x, low Q 2 measurements for Q 2 > 5 GeV 2 show scaling in 
p and a. Thus double asymptotic scaling is a dominant feature of F2 in this region. 



7.4 Extraction of the Gluon Density 

In a QCD analysis the evolution equations were solved numerically in the NLO order ap- 
proximation following the procedure described in [ 53 1 . The splitting functions and the strong 
coupling constant a s (Q 2 ) are defined in the MS factorization and renormalization schemes. 
In the DGLAP evolution equations only three light quark flavours are taken into account. 
Heavy quark contributions are dynamically generated using the photon-gluon fusion (PGF) 



prescription given in [54, |5q1 , extended to NLO according to |pq |. The scale of the PGF 
process has been taken as \/Q 2 + 4m 2 . with a charm quark mass of m c = 1.5 GeV. An un- 
certainty of the charm quark mass of 0.5 GeV was considered which leads to a few percent 
variation of the gluon density. The small contribution of beauty quarks has been neglected. 



The gluon g, the valence quark u v and d v and the non-strange sea S (S = u + d) distri- 

2 




butions are parametrized at Qq = 5 GeV 2 in the following way: 



xg(x) 


= A g x B °(l- 


x) Cg , 




xu v (x) 


= A u x B ^(l- 


x) Cu 


{1 + D u x + E u ^), 


xd v {x) 


= A d x Bd {\- 


x) Cd ( 


1 + D d x + E d ^x~), 


xS(x) 


= A s x Bs (l- 


xf s 


(1 + D s x + E S V^) 



(14) 

The quark and antiquark components of the sea are assumed to be equal, and u is set equal 
to d. As determined in |57|] , the strange quark density is taken to be 5/4. With these 
definitions the proton structure function F2, for nj = 3 and to leading order, is given as 

F2(x, Q 2 ) = 7^x5 + ^-xu v + \xd v . (15) 
to y y 



The normalizations of the valence quark densities are fixed using the counting rules 
Jq 1 u v {x)dx = 2 and Jq 1 d v (x)dx = 1. The normalization A g of the gluon density is obtained 
via the momentum sum rule. Since no isoscalar data are available yet in the small x domain, 
B u = B d is assumed. The parameters B$ and B g which govern the small x behaviour of 
F2 and of the gluon are allowed to be different. For A the value of 263 MeV is taken, as 
determined in [52|. 
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x = 0. 000032 
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=d. ; D002 
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• H1 
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H1 NLO fit 




QVGeV 2 



Figure 11: F 2 (x, Q 2 ) measured by HI together with BCDMS ^Tj and NMC |5fj/ hxed target 
results. The full line corresponds to the NLO QCD fit, see sect. 7. 4, which includes the data 
for Q 2 > 5 GeV 2 . The extension of the curves below 5 GeV 2 represents only the backward 
evolution of the fit. The F 2 values are plotted in a linear scale adding a constant c(x)= 0.6(i- 
0.4) where i is the x bin number starting at i =1 from x = 0.32. The inner error bar is the 
statistical error, the outer corresponds to the full error resulting from adding the statistical 
and systematic error in quadrature. Some HI data points at lower Q 2 where shifted to nearby 
x values for graphical representation of the data. 
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Experiment 


HI 

nvtx 


HI 

svtx 


NMC-p 


NMC-D 


BCDMS-p 


BCDMS-D 


total 


data points 


157 


16 


96 


96 


174 


159 


698 


X 2 (unco, err.) 


174 


13 


157 


153 


222 


208 


931 


X 2 (full error) 


85 


6 


120 


114 


122 


140 


591 


normalization 


1.00 


1.04 


1.00 


1.00 


0.97 


0.97 





Table 5: For each experiment are given: the number of data points used in the QCD fit, the 
X 2 obtained as described in the text using only the uncorrelated errors, the x 2 computed from 
the same fit using the full error on each point and the normalization factors as determined 
from the ht. The HI nominal vertex and shifted vertex data samples are denoted as nvtx 
and svtx respectively. 



In order to constrain the valence quark densities at high x, proton and deuteron results 
from the BCDMS and NMC experiments are also used. To avoid possible higher twist effects, 
data in the ranges Q 2 < 5 GeV 2 , and Q 2 < 15 GeV 2 for x > 0.5 are not included in the fit. 
The small contribution of large rapidity gap events in the HERA data is considered to be part 
of the structure function, as there is no evidence that the QCD evolution of the diffractive 
part of F2 is significantly different from that of the total inclusive F2 . 

The parton densities are derived from a fit of the evolution equations to the data using 
the program MINUIT. For the calculation of the x 2 which was minimized, the statistical 
errors were combined in quadrature with those systematic errors which are uncorrelated. 
For BCDMS only statistical errors were included. In addition a term was added to the 
X 2 to permit variation of the relative normalization of the different data sets. The following 
normalization errors were taken into account: HI (nominal vertex sample): 1.5%, HI (shifted 
vertex sample): 3.9%, BCDMS: 3%, and NMC: 2.5%. The x 2 obtained in this procedure and 
the x 2 computed when considering the full error of each point are given in Table |. 

The result of the fit is shown in Fig. [7| versus x and Fig. 11 versus Q 2 . The fit gives 
a good description of all data used. Only small adjustments of the relative normalizations 
(given in Table ||) are required demonstrating remarkable agreement between these different 
experiments. In Fig. the steep x behaviour of F2 is seen to be described very well by the 
fit. Note that the data for Q 2 < 5 GeV 2 , which were excluded from the fit, are still well 
reproduced by the fit evolved backwards in Q 2 . However, a definite test of perturbative QCD 
in this region requires more accurate data and a study of possible higher twist effects, which 
is beyond the scope of this analysis. The Q 2 dependence at fixed x is also described well over 



the nearly 4 orders of magnitude in Q 2 covered by the HI data, see Fig. 11 . The parameters 
of the initial distributions are listed in Table |6| There are sizeable correlations between these 
parameters which were not studied in detail as the basic aim of this analysis was to extract 
the gluon density. 



Fig. 11 shows the NLO gluon density xg(x,Q 2 ) at Q 2 = 5 GeV 2 and Q 2 = 20 GeV 2 . 
Note that there are no F2 measurements below 5 • 10 at Q 2 = 20 GeV 2 , but in that 
region the gluon is constrained by the data at lower Q 2 via the QCD evolution equations. 
The experimental error band was determined in two steps: the initial error was obtained 
directly from the fit when considering only the uncorrelated errors of the data points which 
are dominated by the statistical errors of the F2 measurement. Then the systematic errors 
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A 9 


B 9 


c 9 






A u 


B u 


C u 


^« 


E u 


2.24 


-0.20 


8.52 






2.84 


0.55 


4.19 


4.42 


-1.40 


A s 


B S 


Cs 


D s 


Eg 


A d 


5d 


Cd 


D d 


E d 


0.27 


-0.19 


1.66 


0.16 


-1.00 


1.05 


0.55 


6.44 


-1.16 


3.87 



Table 6: The values of the parameters at Q 2 =5 GeV 2 of the gluon, the sea quark and 
valence quark densities, as determined from the QCD fit. 

HI 1994 




Figure 12: The gluon density xg(x) at Q 2 = 5 GeV 2 and Q 2 = 20 GeV 2 extracted from a 
NLO QCD fit. The procedure to derive the error bands is explained in the text. 
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introducing point-to-point correlations, as for instance a possible shift in the scattered electron 
energy, were treated in a procedure described in [51]. Their effect was added quadratically 
to the effect from the uncorrelated ones, determining the full experimental error band of the 
measurement of xg. A variation of A by 65 MeV [52] gives a change of 9% on the gluon 
density at 20 GeV 2 which has not been added to the errors shown in Fig. |l^. 

The accuracy of this determination of xg is better by about a factor of two than the HI 
result based on the 1993 data [51]. A rise of the gluon density towards low x is observed 
which is related to the behaviour of F2 oc x~ x . Accordingly, the rise of xg towards low x 
increases with increasing Q 2 . 



8 Summary 

A measurement has been presented of the proton structure function F2(x, Q 2 ) in deep- inelastic 
electron-proton scattering at HERA with data taken in the running period of 1994. The 
integrated luminosity is 2.7 pb _1 which represents a tenfold increase in statistics compared 
to the Fi analysis based on the 1993 data of the HI experiment. The structure function 
measurement includes data from different detector components and running configurations. 
Low Q 2 values are reached using data with the ep interaction vertex shifted from the nominal 
position and with radiative events. The data cover a kinematic range for Q 2 between 1.5 and 
5000 GeV 2 and x between 3.0 • 10~ 5 and 0.32. 

The F2 values presented are obtained using different methods to reconstruct the inclusive 
scattering kinematics. At high values of the scaling variable y > 0.15, due to its superior 
resolution, the method used is based on the scattered electron energy and angle. Lower 
y values are covered with a method which combines electron and hadronic information to 
reduce radiative corrections and calibration errors. A smooth transition is observed from the 
fixed target high x data to the HERA low x data. 

The rise of the structure function with decreasing x at fixed Q 2 is confirmed. The rate 
of growth increases with increasing Q 2 which has been one of the very first predictions of 
perturbative QCD. Approximate scaling in double logarithmic variables depending on x and 
Q 2 is observed using a recent 2 loop QCD calculation. The data are well described in the 
full x and Q 2 range by a NLO fit based on the conventional DGLAP evolution equations. 
The fit results are used to measure the gluon distribution with improved precision down to 
x = 10 -4 . The gluon density rises significantly for decreasing values of x. 
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Table 7: Proton structure function ^(x, Q 2 
The normalization uncertainty, not included 
GeV 2 and 3.9% for Q 2 < 8.5 GeV 2 . 
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) with statistical and systematic errors, part I. 
in the systematic error, is 1.5% for Q 2 > 8.5 
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Table 8: Proton structure function F%{x, Q 2 ) 
The normalization uncertainty, not included i: 
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with statistical and systematic errors, part II. 
i the systematic error, is 1.5%. 
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